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A report is presented on our continued effort to elucidate the continuum limit of Bk using the quenched 
Kogut-Susskind quark action. By adding to our previous simulations one more point at /3 = 6.65 employing a 
56^ X 96 lattice, we now confirm the expected O(a^) behavior of Bk with the Kogut-Susskind action. A simple 
continuum extrapolation quadratic in a leads to _Bk(NDR, 2 GeV) = 0.598(5). As our final value of Bk in the 
continuum we present _B/f (NDR, 2 GeV)=0. 628(42), as obtained by a fit including an Q-jy^(l/a)^ term arising 
from the lattice-continuum matching with the one-loop renormalization. 



1. Introduction 



2. Perturbative matching 



The problem of taking the continuum limit re- 
mains to be the most important issue in the lat- 
tice calculation of Bk using the quenched Kogut- 
Susskind quark action. We have been making ef- 
fort to elucidate the effect of scaling violation over 
two years since 1995. Our early study[|l[ carried 
out at /3 = 5.85 — 6.2 showed a linear decrease 
of Bk in a in contrast to 0{a^) scaling violation 
theoretically predicted The run was then ex- 
tended to P = 6.4[^, at which a departure from a 
linear behavior was observed, indicating the on- 
set of O(a^) behavior. This implied a cancellation 
between an and higher order terms that leads 
to an apparent 0{a) behavior for lower values of 
/3. We have now extended the run to f3— 6.65 em- 
ploying a 56'^ X 96 lattice to settle the issue of the 
continuum extrapolation. We also briefly address 
the problem of 0{al) effect, on which we gained 
insight after the present conference. 



We have slightly revised the method of anal- 
ysis in obtaining Bk in the continuum since 
Lattice 96 Lattice values of Bk are con- 
verted to the continuum value in the MS scheme 
with the naive dimensional regularization (NDR) 
by applying one-loop renormalization|Q,^ at the 
matching scale g* = l/a[Hj^. The one-loop renor- 
malization factor is evaluated with the 3-loop run- 
ning coupling constant ctjjg{q*) with A-^g-g — 0.23 
GeV, estimated in the continuum limit from our 
results for the p meson mass. 

The continuum value of Bk at the physical 
scale p, = 2 GeV is calculated from q* via the 
2-loop running of the continuum renormalization 
group: _ 
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Table 1 



Run parameters and results for i3/(-(NDR, 2GcV) 



p 




La 


BK(NDR,2GeV) 






(fm) 


non-inv. 


inv. 


5.7 


12^24 


2.8 


0.846(1) 


0.822(1) 


5.85 


16^32 


2.3 


0.780(2) 


0.756(2) 


5.93 


20^40 


2.5 


0.752(2) 


0.723(2) 


6.0 


24^64 


2.5 


0.715(2) 


0.683(2) 




18^64 


1.9 


0.717(7) 


0.679(7) 




32^48 


3.1 


0.713(1) 


0.679(2) 


6.2 


32^64 


2.4 


0.662(5) 


0.624(5) 


6.4 


40^96 


2.3 


0.643(7) 


0.607(7) 




32^96 


1.8 


0.658(12) 


0.613(11) 




48^96 


2.7 


0.642(5) 


0.607(5) 


6.65 


56^96 


2.3 


0.635(7) 


0.606(7) 



This matching procedure leaves aj^{l/a)'^ un- 
certainty, which however should vanish logarith- 
mically in the continuum limit 0,|l^. 



3. Run parameters and results 

Table |l| lists the parameters of all our runs to 
date, together with ^/^(NDR, 2 GeV) obtained 
with gauge non-invariant and invariant operators 
at each f3. 

Finite-size studies have been carried out at /3 = 
6.0 and 6.4[P. The finite size effect is less than 
1% above 2.3 fm at both /?, though it increases to 
2% below 2.0 fm at /3 = 6.4. We have therefore 
focused our runs on a physical spatial lattice size 
of La w 2.3 — 2.8 fm. We expect the finite size 
effect to be less than our statistical errors. 

4. Continuum extrapolation 

In Fig. I we present ^/^(NDR, 2 GeV) as a 
function of TOpO. The leftmost points are the new 
data taken at /? = 6.65 on a 56'^ x 96 lattice. 
The new data points strongly support an O(a^) 
behavior, the onset of which was noted in Lattice 
96 at /3 = 6.4. The effect of higher order terms 
is not discernible for /3 > 5.93. Therefore, we 
fit the five points above (3 = 5.93 with the form 
Bk — Co + ci(mpa)^, which is shown by dashed 
lines in Fig. |l|. The continuum extrapolation gives 
Sif(NDR, 2 GeV)=0.616(5) for the gauge non- 
invariant operator, and 0.580(5) for the invariant 
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Figure 1. Gauge non- invariant (circles) and invari- 
ant (diamonds) iJif (NDR, 2GeV) as a function of 
TOpfl, together with a simultaneous fit for the two 
operators including term (solid lines) and sep- 
arate fits quadratic in a(dashed lines) to the five 
pairs of data points for (3 > 5.93. 

one, the average of the two being 0.598(5). 

5. Operator dependence 

It has been noted that the two operators yield 
different values in the continuum[Q (see Fig. |^). 
After Lattice 97 a further analysis was made on 
this point, which we report in this write-up. The 
difference between gauge non-invariant and in- 
variant operators should receive not only O(a^) 
scaling violation but also ctjfgil*)'^ errors from 
the matching procedure. Figure || plots this dif- 
ference as a function of nipa. Small errors result- 
ing from a correlation between the two operators 
allow us to fit the five data points with the form 
bi{j7ipa)'^ + b2aj^{q*)'^ , giving 6i = -0.23(2) and 
62 = 1.73(5) with xVd.o.f = 2.2. The solid 
line indicates the fit, and others show the break- 
down into the a^ (dotted line) and (dashed line) 
contributions. Allowing for a non-zero constant 
60, similar fitting gives bo = —0.032(16), bi = 
-0.44(11) and 62 = 3.4(8). 

Encouraged by this analysis we attempt to fit 
the five points at /3 > 5.93 simultaneously for 
both operators including their correlations, with 
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Figure 2. Difference of (NDR, 2GeV) between 
gauge non-invariant and invariant operators as a 
function of rripa. The solid line represents a fit 
with and c? terms, while the dotted (dashed) 
line is the contribution from the (a^) term. 



the form 



Co + cia 



'C2a-g-g{q*Y and 



Bj^^- = do + dia^ d-2ajjg{q* f . This fit yields 
Co — 0.67(6) and do = 0.71(7), and hence we 
impose the constraint co = do in our final fit. In 
the continuum limit the fit (solid lines in Fig. |^) 
gives Bk (NDR, 2 GeV) = 0.628(42). The error 
is roughly ten times larger than the one from the 
quadratic fit. This originates from uncertainties 
of the coefficient of terms: C2 = —0.5(2.0) and 
d2 = —2.2(2.0). The difference, however, is well 
constrained; C2 — d2 = 1.7 agrees well with 62 — 
1.73 obtained above. 

We find larger coefficients C2 = —1.0(4.2) and 
d2 = —4.3(4.2) when q* — ir/a is used, or 
C2 = 1.6(1.5) and d2 = —3.2(1.5) if mean-field 
improvement of the operators is not made. This 
supports the argument of Ref . ||ll[] . 

The final value depends only weakly on our 
choice of K-jjg = 0.23 GcV: e.g., ^^(NDR, 
2 GeV)=0.627(42) for A^=0.22 GeV and 
0.628(41) for 0.24 GeV in the continuum hmit. 



6. Conclusions 

As our final value of Bk in the continuum limit 
we adopt the result from the fit including the 



term, (NDR, 2GeV) = 0.628±0.042, which in- 
cludes a systematic error due to the 2-loop uncer- 
tainty. The size of the quoted error is 6.6%, which 
roughly equals 3 x ot-jjg{q* = 1/a)^ at our small- 
est lattice spacing 1/a = 4.87GeV at /? = 6.65 
where ap-g(4.87GeV) = 0.147. This magnitude 
of error is unavoidable unless a two-loop calcula- 
tion is carried out for the lattice renormalization. 

Our Bk is consistent with the JLQCD 
value obtained using the Wilson quark action, 
Bif(NDR,2GeV) = 0.562 ± 0.064, in which 
the operator mixing problem is solved non- 
perturbatively with the aid of chiral Ward iden- 
tities |jl|l. 
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